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Using  ^th  a  bare-atom  and  dressed-atom  picture  of  atom-field  interactions,.**^  present  a  theory 
of  the  pressure-induced  resonances  that  can  arise  when  four  radiation  fields  interact  with  a  three- 
level  quantum  system.  The  four  fields  drive  two  coupled  atomic  transitions  which  share  a  common 
level,  the  levels  being  arranged  in  a  “ladder”  scheme.  As  the  field  frequencies  are  varied,  certain 
resonant  structures  in  the  excitation  spectrum  of  the  upper  state  population  appear  only  in  the  pres¬ 
ence  of  collisions.  The  resonant  structures  result  from  processes  that  have  been  referred  to  in  the 
past  as  collisionally  aided  radiative  excitation  (CARE)  and  pressure-induced  extra  resonances  in 
four-wave  mixing  (PIER4).  Wr^liow  tlilf^th  CARE-  and  PIER4-type  resonances  have  a  common 
origin,  which  can  be  linked  to  conservation  of  energy.  In  contrast  to  conventional  PIER4  signals, 
line  shapes  calculated  in  this  work  result  from  a  “single  atom^^rather  than  a  collective  emission. 
Moreover,  the  signals  vanish  unless  the  applied  fields  are  relatively  coherent,  again  in  contrast  to 
conventional  PIER4.  f-' 


1.  INTRODUCTION 

Spectroscopic  line  shapes  provide  us  with  most  of  the 
data  we  have  concerning  atomic  level  spacings  and 
atomic-state  lifetimes.  In  an  atomic  vapor,  the  spectro¬ 
scopic  line  shapes  are  modified  as  a  result  of  collisions 
occurring  within  the  vapor.  The  most  commonly  known 
effect  of  this  type  is  the  so-called  pressure  broadening  of 
spectral  lines,  in  which  absorption  or  emission  lines  are 
broadened  as  a  result  of  collisions.  Perhaps  somewhat 
less  familiar  is  a  whole  series  of  processes  in  which  col¬ 
lisions  result  in  new  structure  in  either  linear  or  nonlinear 
spectroscopic  line  shapes. 

The  collision-induced  features  can  occur  in  reactions  in 
which  both  colliding  atoms  change  their  internal  states  as 
a  result  of  the  combined  collisional-radiative  interaction.’ 
They  can  also  occur  when  a  forbidden  transition  becomes 
allowed  as  a  result  of  collisions.^  In  this  paper,  neither  of 
these  cases  is  discussed.  Instead,  we  consider  the  interac¬ 
tion  of  atoms  with  one  or  more  radiation  fields  that  drive 
allowed  transitions  in  these  “active”  atoms.  The  active 
atoms  undergo  collisions  with  ground-state  perturber 
atoms  which  remain  in  their  ground  state  and  interact 
negligibly  with  the  radiation  fields.  Nevertheless,  the  ac¬ 
tive  atom-perturber  collisions  can  modify  the  active 
atom-radiation  field  interaction  in  a  manner  that  leads  to 
new  features  in  the  absorption  and  emission  line  shapes. 

Collisional  effects  of  this  nature  have  been  broadly  di¬ 
vided  into  two  categories.  The  first  category  can  be 
termed  collisionally  aided  radiative  excitation  (CARE),^ 
and  is  illustrated  in  Fig.  1 .  A  weak  laser  field  having  fre¬ 
quency  n  drives  a  transition  between  the  ground  state  1 
and  excited  state  2  of  an  active  atom.  The  magnitude  of 
the  detuning  A  =  «  — H  («=l-2  transition  frequency)  is 
larger  than  the  decay  rate  of  level  2.  A  second  weak  laser 


has  its  frequency  IT  tuned  in  the  neighborhood  of  the  2-3 
transition  frequency  w'  and  the  population  of  state  3,Piy, 
is  monitored  as  a  function  of  A'=ct)'  — H'.  In  the  absence 
of  collisions,  there  is  a  single  resonance  centered  at 
A'=  — A  w').  When  collisions  occur,  how¬ 

ever,  there  is  an  additional  resonance  centered  at 
A'=0  This  type  of  collision-induced  feature 

can  also  be  observed  in  fluorescence,  where  it  goes  under 
the  name  of  collisional  redistribution  ^  In  that  case,  the 
fluorescence  spectrum  from  level  2  consists  of  a  single 
peak  centered  at  the  laser  frequency  n  in  the  absence  of 
collisions.  With  collisions  present,  there  is  an  additional 
peak  centered  at  the  atomic  transition  frequency  a.  In 
both  cases  above,  the  width  of  the  additional  component 
is  determined  by  the  decay  rate  of  the  optical  coherence 
associated  with  the  transition  levels.  A  physical  explana¬ 
tion  of  these  CARE  resonances  is  given  in  Sec.  II.  At 
this  point,  however,  we  note  that,  since  the  collision- 
induced  features  are  monitored  from  atomic  state  popula¬ 
tions,  they  are  essentially  “single  atom”  effects.  No 
cooperative  emission  is  needed. 

The  second  category  of  collision-induced  features 
occurs  for  the  atom-field  interaction  indicated  in  Fig.  2. 
Two  laser  fields  having  frequencies  H  and  H  -l~  8  drive  the 
1-2  transition.  These  fields  result  in  a  modulation  of  the 
population  of  state  2  at  rate  8.  When  a  third  field  is  ap¬ 
plied  having  frequency  H,  it  is  possible  to  obtain  a  four- 
wave-mixing  signal  which  propagates  with  frequency 
11  —  8.  For  detunings  |A|  =  |(u  — fl|  much  larger  than  the 
natural  width  of  the  1-2  transition  (and  any  Doppler 
widths),  it  is  found  that,  as  8  is  varied,  a  resonant  struc¬ 
ture  centered  at  8  =  0  having  width  [full  width  at  half 
maximum  (FWHM)]  2/2  appears  in  the  signal  field  only 
in  the  presence  of  collisions.  Such  pressure-induced  extra 
resonances  in  four-wave  mixing  (PIER4)  were  predicted 
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FIG.  1.  Energy-level  scheme  for  collisionally  aided  radiative 
excitation  (CARE).  The  detuning  A  =  (u  — ft  is  held  fixed  as  the 
detuning  A'  =  (u'  — ft'  is  varied.  The  upper  state  population  exci¬ 
tation  spectrum  has  a  resonant  structure  centered  at  A'=0  only 
in  the  presence  of  collisions. 

and  observed  by  Bloembergen  and  co-workers.^  ’  They 
differ  from  the  first  category  of  resonances  in  that  their 
width  is  determined  by  an  atomic-state  population  decay 
rate  rather  than  a  coherence  decay  rate.  Moreover,  as 
observed  in  four-wave  mixing,  they  result  from  a 
cooperative  phase-matched  coherent  emission. 

It  is  the  purpose  of  this  paper  to  compare  the  two 
categories  of  resonances  and  to  examine  the  origin  of 
both  categories  of  resonances  using  a  unified  approach. 
As  such,  we  will  show  that  PlER4-type  resonances  can  be 
observed  using  level  schemes  similar  to  that  shown  in 
Fig.  1.  In  other  words,  the  PIER  signals  can  be  observed 
as  a  single  atom  effect  rather  than  in  cooperative  emis¬ 
sion.  Moreover,  the  signals  which  are  calculated  for  the 
level  scheme  of  Fig.  1  will  be  seen  to  depend  on  the  rela¬ 
tive  coherence  of  the  applied  laser  fields,  in  contrast  to 
the  conventional  PIER4  signals. 

This  paper  is  organized  as  follows.  In  Sec.  II  a  physi¬ 
cal  description  of  both  CARE  and  PIER4  is  given  and 
similarities  and  differences  between  the  two  categories  of 
resonances  are  noted.  In  Sec.  Ill  we  use  a  three-level 
atom  and  four  laser  fields  to  derive  a  signal  that  possesses 


^  1 
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FIG.  2.  A  four-wave-mixing  signal  having  wave  propagation 
vector  k4  =  k| -t-k,  — k,  and  frequency  ft4  =  ft  — 8  is  produced 
when  three  fields  having  propagation  vectors  ki.kj.ki  and  fre¬ 
quencies  ft|  =  ft,  ft;  =  ft +  8,  ft,  =  ft  drive  a  two-level  atomic 
transition.  A  resonant  structure  in  the  signal  intensity  centered 
at  8  =  0  appears  only  in  the  presence  of  collisions  (pressure- 
induced  extra  resonance  in  four  wa\c  rr.i>iiiL’  ;riCR4)j. 


both  CARE-  and  PIER4-type  features.  The  common  ori¬ 
gin  of  both  resonant  structures  emerges  in  this  section. 
The  calculation  in  Sec.  Ill  is  carried  out  using  a  semiclas- 
sical  approach.  In  Sec.  IV  this  calculation  is  repeated  us¬ 
ing  a  fully  quantized,  dressed-atom  approach.  It  is  ar¬ 
gued  that  the  pressure-induced  resonances  can  be  ex¬ 
plained  on  the  basis  of  conservation  of  energy.  In 
planned  companion  papers,  we  shall  investigate 
pressure-induced  signals  which  arise  in  (a)  the  fluores¬ 
cence  beats  emitted  by  a  two-level  atom  driven  by  two 
laser  fields;  (b)  the  excitation  of  a  “four-level”  atom  using 
four  applied  laser  fields;  and  (c)  the  two-photon  ionization 
of  an  atom  produced  using  two  laser  fields.  In  those  pa¬ 
pers  a  semiclassical  dressed-atom  picture  will  be  intro¬ 
duced  which  provides  an  alternative  approach  for  carry¬ 
ing  out  the  calculations  than  the  methods  developed 
herein. 

II.  DISCUSSION  OF  PRESSURE-INDUCED  EXTRA 
RESONANCES  (PIER) 

In  this  section  we  show  that  the  pressure-induced  extra 
resonances  observed  in  collisionally  aided  radiative  exci¬ 
tation  (CARE)  and  four-wave  mixing  (PIER4)  are  inter¬ 
related.  The  existence  of  both  of  these  pressure-induced 
extra  resonances  can  be  traced  to  a  collisional  interaction 
in  which  part  of  the  internal  energy  of  the  (active-atom 
plus  radiation  field)  is  transformed  into  translational  en¬ 
ergy  of  the  colliding  atoms.  Both  a  standard  “bare- 
atom”  picture  and  a  dressed-atom  picture  is  used  to  illus¬ 
trate  this  aspect  of  the  theory.  It  is  always  assumed  that 
level  1  is  the  ground  state  and  that  collisions  alone  are 
insufficiently  energetic  to  cause  transitions  between  the 
internal  atomic  states  under  consideration. 

A.  CARE 

The  level  scheme  for  a  typical  CARE  reaction  is  shown 
in  Fig.  1.  One  laser  has  its  frequency  ft  held  fixed  at  a 
value  which  gives  a  detuning  |A|  =  |w  — ftl  that  is  large 
compared  with  the  natural  and  Doppler  widths  associat¬ 
ed  with  the  1-2  transition,  but  small  or  comparable  to  the 
inverse  duration  of  a  collision.  The  Rabi  frequency  asso¬ 
ciated  with  this  field,  which  drives  the  1-2  transition  only, 
is  denoted  by  x  (;(f=/ 12^/2^;/ 12  is  the  dipole  moment 
matrix  element,  E  is  the  field  amplitude).  A  second  laser 
field  (frequency  ft';  Rabi  frequency  x'  /2fi)  drives 

the  2-3  transition  only.  As  the  detuning  A'  =  w'  — ft'  is 
varied,  the  population  of  state  3  is  monitored.  In  this,  as 
well  as  in  all  subsequent  discussions,  it  is  assumed  that  the 
laser  fields  are  weak;  results  are  calculated  to  lowest  non¬ 
vanishing  order  in  the  field  amplitudes. 

A  question  that  may  be  posed  is  for  what  detunings  A' 
does  resonant  excitation  of  level  3  occur?  The  answer  to 
this  question  at  first  appears  to  be  paradoxical.  Consider 
first  a  density-matrix  treatment  of  the  problem.  Two 
types  of  density-matrix  chains  lead  from  an  initial  popu¬ 
lation  Pi  I  to  a  final  population  p,,  (p,^  is  the  i-j  density- 
matrix  element).  The  first  of  these  involves  the  inter¬ 
mediate  state  population  p,,  and  can  be  represented 
.schemaiicaily  as 
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Pii  P22  P33 

/i*  /r* 

P21  Pn 


(1) 


The  Rabi  frequencies  in  chain  (1)  specify  the  field  that 
produces  a  given  density-matrix  element  from  a  preced¬ 
ing  one  in  a  perturbative  solution  to  the  problem.  For 
chain  (1),  which  can  be  referred  to  as  a  stepwise  (SW) 
chain,  the  population  P33  varies  as  Moreover, 

the  excitation  appears  as  the  sequential  absorption  of  a 
photon  from  field  1  and  a  photon  from  field  2.  Conse¬ 
quently,  there  should  be  a  resonance  centered  at 
A'  =  <u'  — n'=0  for  the  second  step  in  the  absorptive  pro¬ 
cess. 

The  second  perturbative  chain  leading  to  a  contribu¬ 
tion  to  P33  of  order  is  represented  schematically 


by  the  chain 

r* 

X 

... 

\ 

Pu 

P33  ■ 

(2) 

X, 

/r- 

P21 - p,,— 

V 

- P32 

V* 

This  chain  can  be  termed  a  two-quantum  (TQ)  chain 
since  it  requires  the  simultaneous  action  of  both  fields. 
The  presence  of  p,3  (or  P3,)  in  the  chain  indicates  a  two- 
quantum  resonance  condition  (or 

A-l-A'=0).  Thus,  from  a  density-matrix  approach,  there 
appears  to  be  two  resonances,  a  A'=0  resonance  from 
chain  (1)  and  a  A'  =  —  A  resonance  from  chain  (2). 

On  the  other  hand,  consider  what  happens  if  the  prob¬ 
ability  amplitude  of  state  3  is  calculated.  In  an  amplitude 
picture,  state  2  is  coupled  to  state  1  by  field  X  state  3 
to  state  2  by  field  x'-  When  the  amplitude  calculation  is 
carried  out  (and  all  transients  have  died  away),  only  the 
two-quantum  resonance  ( A'  =  —  A )  remains.  There  seems 
to  be  a  paradox — one  resonance  from  an  amplitude  ap¬ 
proach  and  two  from  a  density-matrix  approach. 

The  resolution  of  this  paradox  lies  in  the  limited  appli¬ 
cability  of  the  amplitude  approach.  The  amplitude  argu¬ 
ment  given  above  is  valid  only  in  the  absence  of  col¬ 
lisions.  If  collisions  occur,  the  state  amplitudes  undergo 
rapid  phase  changes  at  the  time  of  the  collision,  leading 
to  new  contributions  for  the  final-state  population.  It  is 
easier  to  incorporate  the  collisional  effects  into  the  densi¬ 
ty  matrix  rather  than  the  amplitude  equations.  Collisions 
result  in  a  decay  of  off-diagonal  density-matrix  elements, 
with  no  population  decay.  In  the  amplitude  picture,  one 
must  consider  the  role  of  a  single  collision  at  each  active 
atom  site,  sum  over  the  contributions  to  the  final-state 
density  matrix  element  from  each  site,  and  then  average 
over  collision  histories. 

The  question  regarding  the  number  of  resonances  can 
now  be  clarified.  In  the  absence  of  collisions,  the  ampli¬ 
tude  approach  is  valid  and  there  is  a  single  resonance 
centered  at  A' =— A  (  O -co '  —  ,4- This  result  can  be 
further  understood  by  making  reference  to  Fig.  1.  The 
fields  are  assumed  to  be  monochromatic,  and  absorption 


of  the  fields  occurs  by  atoms  with  a  well-defined  ground- 
state  energy.  Imagine  for  the  moment  that  state  3  was 
metastable.  Then  owing  to  conservation  of  energy,  exci¬ 
tation  to  level  3  could  occur  only  when  fl  +  fl'=w  +  oi\ 
independent  of  the  width  of  level  2.  When  level  3  has  a 
finite  width,  the  conclusion  is  not  altered  except  that  the 
resonance  condition  need  be  satisfied  only  within  the  nat¬ 
ural  width  of  level  3.  The  existence  of  a  single  resonance 
is  linked  to  conservation  of  energy.^  How  is  this  result  ex¬ 
plained  in  the  density-matrix  approach?  It  turns  out 
that,  in  the  absence  of  collisions,  the  TQ  chain  (2)  con¬ 
tains  a  term  which  exactly  cancels  the  A'=0  resonance  of 
the  SW  chain  (1).*  This  result  can  be  viewed  as  a  destruc¬ 
tive  interference  effect  between  the  two  chains. 

When  collisions  are  present,  the  amplitude  approach 
discussed  above  is  no  longer  valid.  Moreover,  is  not  pos¬ 
sible  to  use  a  conservation-of-energy  argument  to  rule  out 
the  possibility  of  a  A'  =  0  (n'=w')  centered  resonance, 
since  the  colliding  atoms  can  now  compensate  for  the  de¬ 
tuning  A  by  a  corresponding  change  of  fiA  in  their 
translational  energy.  In  other  words,  state  2  can  be  popu¬ 
lated  by  a  CARE  reaction  in  which  the  energy  mismatch 
^A  is  provided  by  this  collision.  Absorption  from  level  2 
then  results  in  the  SW  resonance  centered  at  Cl'  =  co'.  As 
a  consequence  of  collisions,  the  contribution  of  the  SW 
chain  grows  linearly  with  the  pressure,  and  the  Il'  =  w' 
centered  resonance  is  no  longer  canceled  by  the  TQ 
chain.  With  collisions  present,  both  the  A'=0  and 
A'  =  —  A  resonances  exist. 

The  appearance  of  the  A'=0  resonance  can  be  viewed 
as  the  “destruction  of  destructive  interference.”  The  des¬ 
tructive  interference  of  the  SW  and  TQ  chains,  which 
combined  to  eliminate  the  A'  =  0  resonance  of  the  SW 
contribution  in  the  absence  of  collisions,  is  itself  “des¬ 
troyed”  in  the  presence  of  collisions.  Although  such  an 
interpretation  is  possible,  it  seems  to  us  that  the  argu¬ 
ments  relating  to  energy  tonsci  v<ition  offer  a  inoic  physi¬ 
cal  interpretation  of  the  results. 

Up  to  now,  we  have  employed  what  can  be  termed  a 
bare-atom  picture  (BAP)  of  the  excitation  process.  It  is 
“bare”  in  the  sense  that  eigenstates  of  the  isolated  atom 
are  used.  It  will  prove  useful  to  also  consider  a  dressed- 
atom  picture  (DAP)  of  the  reaction.’' In  the  DAP,  j 
eigenstates  of  the  atom  plus  the  field  driving  the  1-2  tran-' 
sition  serve  as  the  basis  states  for  the  problem.  In  the; 
weak-field  limit,  the  appropriate  eigenstates  and  eigenen-" 
ergies  (£,  is  arbitrarily  taken  equal  to  zero)  are 


I)  =  |l;n,n2>  +  “|2;n,-l,n2>, 


=  0 


(3a) 


Ill)  =  -^|l;ni,n2)  +  l2:n,-l,n2).  ^11 


(3b) 


|III>  =  |3;n,  —  l.Oj “  1 ).  £in=^(A  +  A' 


(3c) 


where  n,  is  the  number  of  photons  in  the  first  field  and  /I2 
is  the  number  in  the  second.  These  dressed  states  are  il¬ 
lustrated  in  Fig.  3.  States  III)  and  llll)  have  energies 
^A  and  ^(A  +  A'),  respectively,  relative  to  state  ll).  As 
A'  is  varied,  state  llll)  moves  vertically  in  Fig.  3,  and 
resonances  can  occur  whenever  state  |I1I>  becomes  de¬ 
generate  with  either  states  |I )  or  III ). 

In  the  absence  of  collisions,  assuming  that  the  field  V  >s 
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FIG.  3.  Dressed-state  energy  levels  in  the  weak-field  limit. 
As  A'  is  varied,  the  energy  of  state  I  III )  moves  vertically  in  the 
figure.  Resonances  can  occur  when  £|„  =  £|  (centered  at 
A'=  —  A)  or  when  £,ji  =£ij  (centered  at  A' =  0).  In  the  absence 
of  collisions,  state  |II)  is  not  populated  and  only  the  A'=— A 
resonance  exists. 


turned  on  in  a  time  long  compared  with  |A|“',  only 
dressed  state  |l)  is  populated.  The  absence  of  state  |II) 
population  is  directly  linked  to  the  conservation  of 
energy — there  is  no  mechanism  to  provide  the  energy  ^A 
to  excite  state  |II).  As  a  consequence,  the  only  reso¬ 
nance  that  occurs  in  the  absence  of  collisions  is  that  for 
which  state  |III)  is  degenerate  with  state  |I),  i.e.,  when 
(A  +  A')  =  0.  When  collisions  occur,  the  situation 
changes.  Collisions,  which  do  not  couple  the  bare  states, 
do  couple  dressed  states  |l)  and  lll>  since  each  of  these 
states  is  a  superposition  of  states  1 1  >  and  1 2 ) .  Asa  result 
of  collisions,  state  !ll)  is  populated  and  a  new  resonance 
centered  at  A'=0  appears  when  state  |III)  is  degenerate 
with  state  |II).  From  either  the  DAP  or  BAP,  the  same 
conclusions  are  reached  concerning  the  existence  of  the 
resonances.  In  both  descriptions,  the  absence  of  the 
A'  =  0  resonance  can  be  predicted  on  the  basis  of  conser¬ 
vation  of  energy  arguments. 

B.  PIER4 

Pressure-induced  extra  resonances  have  been  observed 
using  the  four-wave-mixing  geometry  depicted  schemati¬ 
cally  in  Fig.  2.  Instead  of  considering  the  level  scheme  of 
Fig.  2,  we  choose  to  analyze  that  of  Fig.  4,  since  it  is  simi¬ 
lar  to  that  discussed  above  for  CARE.  Four  fields,  hav¬ 
ing  Rabi  frequencies  lx<X  kX'^X'  and  corresponding 
frequencies  (fl,n  +  5,n',fl'  — 8')  are  incident  on  the 
“three-level"  atom  shown  in  Fig.  4.  Fields  x  and  V  f 
drive  the  1-2  transition  only  and  fields  x'  and  x  drive 
the  2-3  transition  only.  The  detuning  |Al  is  much  larger 
than  the  natural  width  of  level  2  (and  also  larger  than  the 
Doppler  width  asso^i-tcd  with  the  1-2  transition  frequen¬ 
cy),  while  the  detuning  |8l  is  compaiable  with  the  decay 
rate  f  2  2.  In  the  general  case  to  be  considered  in 

Sec.  Ill,  both  8  and  8'  can  be  varied.  For  the  present  dis¬ 
cussion,  however,  8'  is  set  equal  to  8. 

The  frequencies  are  chosen  to  satisfy  an  "energy- 
conservation”  condition  that  is  analcg.'«,.s  *•'  ih:*t  <.>1  <  nn- 
ventional  four-wave  mixing.  In  the  PIER4  experiments, 
three  fields  are  incident  having  frequencies  Il|  =  fl,  fl2 
=  n-)-8,  11,  =  fl.  The  generated  field  has  a  frequency  114 
whose  frequency  is  derived  from  the  energy-conservation 
condition 


FIG.  4.  Energy-level  diagram  and  field  configuration  for  ob¬ 
serving  pressure-induced  extra  resonances  via  a  four-wave- 
mixing  interaction.  Four  fields  arc  incident  on  the  ‘th'-cc-lever’ 
atom  as  shown.  Fields  E  and  E,  drive  the  1-2  transition  and 
fields  E'  and  E'  drive  the  2-3  transition.  As  the  population  of 
level  3  is  monitored  for  fixed  and  11'.  a  resonant  structure  cen¬ 
tered  at  6  =  0  appears  in  the  excitation  spectrum  only  when  col¬ 
lisions  are  present. 

n4=n,-fn,-n2=n-s .  t4) 

[The  fields  must  also  satisfy  momentum  conservation 
k4  =  k|-l-kj  — k2  (k,  is  the  propagation  vector  for  field  1) 
and  phase  matching,  ^4  =  04/0.]  If,  in  Fig.  4,  we  associ¬ 
ate  fields_l,  2,  and  ^with  the  fields  having  frequencies 
n,  =  n,  [l2—il  +  b,  Il,  =  n',  respectively,  then  the  fourth 
field  114  =  H'  —  8  satisfies 

04=12,  + n,--n,=n'-8 ,  (5) 

in  analogy  ..  uh  '  o.  (4).  Thus,  our  fields  satisfy  an  equa¬ 
tion  similar  ,  energy-conservation  condition  of  four- 
wave  mixing,  .v.-,  will  be  seen  below,  however,  the  fact 
that  the  fields  are  chosen  to  satisfy  an  cnergy- 
con.servation  condition  does  not  invalidate  the  conclusion 
that  pressure-induced  extra  resonances  can  be  explained 
in  terms  of  conservation  of  energy  arguments  as  was  done 
above  for  CARE. 

The  discussion  parallels  that  given  for  CARE.  We  first 
adopt  a  bare-atom  approach  (BAP).  In  this  section  wc 
consider  only  the  contribution  to  p,,  (Fig.  4)  pioportional 
to  the  product  of  the  four  Rabi  frequencies  of  the  fields. 
1  he  two  perturbation  chains  which  give  rise  to  this  con¬ 
tribution  are 


and 
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As  in  CARE,  the  SW  chain  (6)  leads  to  resonant  struc¬ 
tures  in  P33  whenever  ft'  [or  (ft'  — 8)]  is  equal  to  w'." 
Moreover,  there  is  a  new  resonant  structure  that  appears 
centered  at  8  =  0  having  width  27^2-  This  resonance  can 
be  viewed  as  arising  from  an  interference  between  the  two 
pathways  leading  to  in  chain  (6).  The  TQ  chain  con¬ 
tains  the  two-photon  resonance  centered  at  ft -I- ft' 
=(o+(o',  having  width  73.  In  the  absence  of  collisions, 
however,  the  resonant  structures  at  ft'  =  (u',  (ft'  — 8)  =  w', 
and  8=0  of  the  SW  chain  are  exactly  canceled  by  the  TQ 
chain,  just  as  in  CARE. 

This  result  can  again  be  explained  using  an  amplitude 
rather  than  density-matrix  approach.  Using  probability 
amplitudes,  one  calculates  the  upper  state  amplitude  as 
the  sum  of  four  terms,  each  representing  a  two-photon 
transition  to  state  3  which  varies  as  the  product  of  a 
lower  (Jt  or  7  +  )  and  upper  (7'  or  x'-  )  field  Rabi  frequen¬ 
cy.  Since  each  term  contains  two-photon  resonance  con¬ 
ditions  only,  when  the  amplitude  is  squared,  the  resultant 
upper  state  population  P33  exhibits  resonances  only  when 
(ft,ft-(-8)  +  (ft',ft'  — 8)  =  <a  +  w'.  This  amplitude  argu¬ 
ment  is  valid  only  if  collisions  are  absent.  Once  collisions 
occur,  this  amplitude  argument  fails  and  one  finds  that 
the  SW  chain  resonances  at  ft'=w',  ft'— 8=«',  and  8=0 
are  not  totally  canceled  by  the  TQ  chain  (7).  Since  the 
vanishing  of  the  SW  resonances  is  implicitly  linked  to  the 
absence  of  a  collision  channel  in  which  internal 
(atom  -1-  field)  energy  can  be  transformed  into  translation¬ 
al  energy  of  the  colliding  atoms,  we  again  see  that  the  ab¬ 
sence  of  PIER  can  be  associated  with  conservation  of  en¬ 
ergy.  This  conclusion  is  strengthened  if  the  reaction  is 
viewed  using  a  dressed-atom  approach  (DAP). 

In  the  DAP  there  are  now  seven  relevant  atom-field 
states, 

|l;/j,,n2:”3>”4).  |2;n,-l,n2:«3."4)  . 

|2;n|,n2“l;«3.”4K  |3;n,  -  1,«2:«3  “  1>”4 )  . 

|3;n  ,,112  -  l;n3  -  l,n4  >  , 

|3;n,,n2“l;«3>”4~l)  - 


•Pl3- 


■P23 


X 


X' 


P33  +  C.c. 


-P31 


-P32 


(7) 


|R>  =  -^|l;«,,n2;«3,n4)  +  l2:ni  -  l,n2;«3.'»4) 

; 

(8b) 

x\ 

|C>-  |l;n,,n2;M3,n4 

)-l-l2;n,,n2-l;n3,n4 

> 

Ec=mA-8)  , 

(8c) 

along  with  the  remaining  eigenkets  and  eigenenergies 

l4>  =  l3:n,  -l,n2;'»3-l.”4^. 

£4=^(A-I-A') 

(8d) 

|5>  =  |3;n,-l,n2;ti3,n4-l). 

£5=^(4 -HA' -1-8) 

(8e) 

|6)  =  l3;n,,n2-l;«3-l,M4). 

(80 

|7>  =  |3;n,,n2“l:'»3.«4~J^. 

E-,=fi{A  +  M  . 

(8g) 

For  the  moment,  we  consider  transitions  to  final  states 
|4)  and  |7)  only  (a  method  for  isolating  these  terms  is 
discussed  in  Secs.  Ill  and  IV).  The  appropriate  energy 
level  diagram  is  shown  in  Fig.  5. 

In  the  absence  of  collisions,  owing  to  conservation  of 
energy,  dressed  states  \B)  and  |C)  are  unoccupied,  so 
that  transitions  to  states  |4)  and  |7)  must  originate  in 
state  \  A).  As  8  is  varied  the  energies  of  states  \  A  ),  l4), 
and  |7)  remain  unchanged.  Consequently,  there  can  be 
no  resonant  structure  centered  at  8  =  0.  Note  that  this  re¬ 
sult  is  now  directly  related  to  conservation  of  energy, 
since  the  absence  of  population  in  states  \B  )  and  jC )  is  a 
consequence  of  energy  conservation. 

When  collisions  occur,  states  i5)  and  lO  are  popu¬ 
lated  and  a  coherence  is  created  between  these  states.  We 
must  now  consider  the  combined  effect  of  the  .4— >-4,7, 
R— *4,7,  and  C— <-4,7  transitions.  As  8  is  varied,  a  reso¬ 
nant  structure  in  the  P33  excitation  spectrum  is  observed 


c 

B 


Ti6 

1 

A-i 


— f - 4,7 

•h(A+A') 


where  the  photon  numbers  refer  to  fields  T,T  +  ,T'.T'-  re¬ 
spectively.  The  fields  x  ^nrl  7+  and  considered  as  the 
“dressing”  fields.  In  the  weak-field  limit,  the  cigenkets 
and  eigenenergies  of  the  (atom -f- dressing  fields)  system 
are  given  by 

\A)  =  \V,nt,n2-,nj,n^)  +  ^\2;ni-l,n2-,n2,n^) 

(8a) 

+  -^^\2■,n^,n2-\■,n2,n^),  E 4=0  , 


FIG.  5.  Energy-level  diagram  for  some  of  the  dressed  states 
associated  with  the  atom-field  interaction  of  Fig.  4.  Fields  E 
and  E,  serve  as  the  dressing  fields.  States  |4>  and  |7)  (which 
are  degenerate  when  6'  =  6)  correspond  to  states  in  which  the 
atom  is  in  level  3.  For  fixed  A  and  A',  excitation  to  states  l4) 
aiid  !?)  o>.cuis  ill  the  wings  of  the  A  *4,7,  B  *4,t,  and 
C  — *4, 7  transitions,  as  shown.  In  the  absence  of  collisions,  only 
state  (  4  )  is  populated.  With  collisions  present,  states  |B  )  and 
|C)  are  populated  and  a  B-C  coherence  is  produced.  As  6  is 
varied,  a  “level-crossing"  resonance  centered  at  8  =  0  appears  in 
the  excitation  spectrum  of  states  |4>  and  |7>,  resulting  from  the 
collision-induced  creation  of  the  B-C  coherence. 
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centered  at  8=0  having  width  2y2-  The  8=0  resonance 
can  be  viewed  as  a  level  crossing  of  states  |B )  and  jC)  in 
the  DAP.'°  The  intensity  of  this  level-crossing  resonance 
is  proportional  to  the  magnitude  of  the  coherence  pg^, 
which  vanishes  (along  with  pgg  and  pcc)  in  tiie  absence 
of  collisions. 

To  conclude  this  section,  we  comment  briefly  on  con¬ 
ventional  PIER4  for  two-level  atoms.  In  PIER4,  energy 
conservation  is  seemingly  built  in  since  condition  (4)  is 
satisfied.  Two  “pump”  photons,  each  having  frequency 
n  are  transformed  into  a  probe  and  signal  photon  having 
frequencies  (n  +  8)  and  (11  —  8),  respectively.  However, 
as  we  have  seen  above,  although  energy  is  conserved  for 
the  fields’  photons,  it  does  not  rule  out  an  explanation 
based  on  energy-conservation  arguments.  If  a  DAP  is 
used,  the  dressed  states  |  /4  ),  |JS  ),  and  |C )  are  still  as  in¬ 
dicated  in  Fig.  5.  What  we  plan  to  show  in  a  future  work 
is  that  the  PIER4  signal  can  be  interpreted  as  a  stimulat¬ 
ed  Raman  transition  between  states  iR)  and  |C)  with 
the  third  field  (at  frequency  fl)  providing  the  Raman 
pump  to  produce  a  Stokes  field  at  frequency  (11  —  8).  The 
physical  arguments  relating  to  the  origin  of  the  PIER4 
resonances  in  such  a  stimulated  Raman  scattering  would 
then  be  analogous  to  those  presented  above  for  three- 
level  atoms  interacting  with  four  fields. 


III.  CALCULATION  OF  PIER  SIGNAL; 
BARE-ATOM  APPROACH 


An  explicit  calculation  of  the  PIER  signal  for  the  level 
scheme  of  Fig.  4  is  now  given.  The  calculation  is  first 
carried  out  using  a  semiclassical  theory  (classical  fields, 
quantum-mechanical  atoms)  and  is  then  repeated  using  a 
fully  quantized  approach. 

Four  fields  are  incident  on  atoms  in  an  atomic  vapor 
having  the  energy-level  structure  shown  in  Fig.  4.  The 
first  two  fields,  which  are  nearly  resonant  with  the  1-2 
transition,  have  electric  field  vectors  given  by 


E(R,t)  =  |€£e'“‘'^-“'’-l-c.c.  , 


2' 

IP 


c  20  I''  c  'Ik+  R-(n  +  8l/] 

E+(R,t)  =  y6+E^+e  c.c.  , 


(9a) 

(9b) 


where  e,e+  are  (complex)  unit  polarization  vectors, 
E,E+  are  (complex)  field  amplitudes,  k,k+  are  field  prop¬ 
agation  vectors,  n,(ft-f8)  are  the  field  frequencies,  and 
c.c.  stands  for  “complex  conjugate.”  The  second  two 
fields,  which  are  nearly  resonant  with  the  2-3  transition, 
have  electric  field  vectors  given  by 

_ I 


E'(R,t)  =  le'£V'"‘''^'““'''-l-c.c.  , 


E'_(R,r)  =  ie'_£:'_e 


,  i[k'_  R-dl'-S'K] 


-f-c.c. 


(9c) 

(9d) 


The  difference  between  the  1-2  and  2-3  transition  frequen¬ 
cies,  denoted  by  a)  and  co',  respectively,  is  sufficiently 
large  to  insure  fhat  fields  E  and  E+  drive  the  1-2  transi¬ 
tion  only  and  fields  E'  and  E'_  drive  the  2-3  transition 
only. 

The  calculation  is  carried  out  in  terms  of  atomic 
density-matrix  elements  p,y  whose  time  development  is 
given  by 

^  ^relaxation  ’  (10) 

where 

H=Ho+y  (11) 

is  the  Hamiltonian  for  the  system  without  relaxation. 
The  quantity  Hq  is  the  free-atom  Hamiltonian  and  V  is 
the  atom-field  interaction  Hamiltonian  taken  to  be  of  the 
form 

F=-/4-E7-(R,t)  ,  (12) 

where /4  is  the  atomic  dipole  moment  operator  and  Ey^  is 
the  total  electric  field,  that  is,  a  sum  of  the  four  fields 
given  in  Eq.  (9).  The  time  derivative  in  Eq.  (10)  is  to  be 
interpreted  as  a  total  time  derivative  equal  to  d/dt  -f  v- V. 

To  account  for  relaxation,  we  adopt  a  rather  simple 
(but  adequate)  model  to  incorporate  the  effects  of  spon¬ 
taneous  emission  and  collisions.  As  a  result  of  spontane¬ 
ous  emission,  level  2  decays  to  level  1  at  rate  Yi  and  level 
3  to  level  2  at  rate  ^3.  Collisions  with  ground-state  buffer 
gas  atoms  are  assumed  to  affect  active-atom  level  coher¬ 
ences  only,  leading  to  a  decay  rate  r,y  =  ry,  for  the 
density-matrix  element  (/#;').  The  rate  F,,  is  propor¬ 
tional  to  the  buffer  gas  pressure.  Recall  that  level  1  is  the 
ground  state,  )'|=0. 

When  the  relaxation  terms  are  included  in  Eq.  (10),  and 
an  interaction  representation  defined  by 


n  =n  p-''kR-n'l 

P12  Pl2^  > 

(13a) 

n  =n  « “'Ik'  R-n’O 

P23  P23^  > 

(13b) 

„  . -iiic+k')R-(n  +  n'i(] 

Pl3“Pl3^  > 

(13c) 

Pij  =p]i 

(13d) 

is  introduced,  one  finds  density-matrix 
evolve  as 

elements  which 

3^22 +  V  •  Vp22  =  -  y  2P22  +  r  3P33  - '  (  T*P2  I  -  TP  12 )  “ '  <  T*  e 


—  i(»rR  — 6f)^  _ V-  x,<(»rR-6/ 

P21  X  +  ^ 


’p.2) 


-'(T'p23-T'*P32)-'(T-e'''''"^*''’p23-T-e 


P32'  ' 


(14a) 


dp33/dt+\Vp,y 

8p,2/3t  -l-v-Vpt,= 


-[(yi2+ri2/-  iiA+k-v)]p,2-i-iAr*(p22-Pii)+'TV 


;(it  R  -  8r 


(P22-Pl 


9P23/9f  +  V-Vp23=-[(y23  +  r23)-t(A'-(-k'-V)]p2j  +  lY*(p33-p22)  +  /y'*C 


-i(A'-R  +  6'f ), 


(P33"‘P22' 


(14b) 

(14c) 


+  'TPi3  +  'T+e'‘''^  ®'’p,3  , 


(14d) 
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ap,3/ar+vVp,3=-((r,3+r,3)-i[(A+A')+(k+k')-v]|p„+i;t*p23+'ne'''''*"®'’P23 


where 

ro=(r,+r;)/2 .  (i5) 

.Y=(2|/|l>-?f:/2«,  Y+  =  <2|/|l>-t  +  £+/2i5  ,  {16a) 

X'={3y.2)  e'E'/2fi,  y'  =<3l/tl2>-6_£'_/2^  (16b) 

A=&)  — fl,  A3.  =6)  — (D  +  6)  (17a) 

A'  =  cj'-n',  A'_=w'-(a'-8')  (17b) 

(c  =  k3.— k,  <c'  =  k'_— k'.  (18) 

In  writing  Eqs.  (14),  it  is  assumed  implicitly  that  the  mag¬ 
nitudes  of  all  the  frequencies  defined  in  Eqs.  (15)- (17)  are 
smaller  than  the  inverse  duration  of  a  collision  (impact 
approximation).  This  assumption  is  needed  to  enable  us 
to  model  collisional  effects  using  simple  decay  parame¬ 
ters.  The  lower  state  population  does  not  appear  explicit¬ 
ly  on  the  left-hand  side  of  Eqs.  (14),  but  is  determined  by 
the  condition 

Pll+^22+/’33=^('')  .  (19) 

where  H^(v)  is  the  equilibrium  velocity  distribution, 

lf'(v)  =  (7ru^)'^^‘e-‘''“' ,  (20) 

in  which  u  is  the  most  probable  atomic  speed. 

We  seek  a  steady-state  solution  of  Eqs.  (14)  for  P33  to 
lowest  nonvanishing  order  in  the  field  amplitudes,  corre¬ 
sponding  to  terms  of  order  lYl%'I^TT*+ lT'I^TT+T'T-. 
etc.  In  this  perturbation  limit,  it  is  sufficient  to  consider 
the  following  truncated  form  of  Eqs.  (14): 

ap22/9f+v-Vp22 

= -r2P2:-'[(T*  +  T+e  '^“®'’)p2i 


-(Y+Y+^''‘*''“®'’)pi2] 


ap33/8t-|-vVp33 


=  -r3P33  +  '[(T'+T-e'*'''’‘’^*''’)p23 

-(Y'*  +  Y-e“''''''  ‘'‘"'’''’)P32]  .  (21b) 

ap,2/at-fvVp,2=-(r'i2-'^ipi2 

-/(Y*+Y  +  e^''''*“®'’)l^(v)  ,  (21c) 

ap23/a/+vVp23 

+  '(T+Y  +  e''*'^''*'')p,3  ,  (21d) 


ap,3/af-)-vVp,3=-[y',3-((A-f-A')]pi 


-i(X'*+X'*e 


-/(ic'  R  +  S'Oi;^ 


)pi2  >  (21e) 


y'll  T>'2‘*‘I'l2>  y23“>'23+ r'23>  +  (22) 

A  =  A-l-k-v,  A'  =  A-fk'-v.  (23a) 

For  future  reference,  we  also  define 

A  3.  =A+ -l-k^-v,  A =A'_ +k'_-v  (23b) 

S  =  6-#c-v,  S'  =  6'-f-»f'-v  ,  (24) 

where  A+  and  A'_  are  defined  in  Eqs.  (17).  In  writing 
Eqs.  (21),  we  set  p|,=  lF(v);  higher-order  corrections  to 
Pi,  lead  to  contributions  to  P33  which  are  down  by  order 
IyP  from  those  calculated  in  this  work. 

To  first  order  in  the  fields,  it  follows  from  Eq.  (21c) 
that  P12  can  be  written  as 


Pl2  =  Pl2(^)+Pl2(^+le 


—  — 6f ) 


When  Eq.  (25)  is  substituted  into  Eq.  (21c),  one  finds 
steady-state  values 

Pi2(A)= -/y*H'(v)/{y',2“'^1  >  (26a) 

P,2(A3, )= -/Y+ l^'(v)/(r'i2“'^+ )  •  (26b) 

Having  obtained  pi2  (along  with  p2i~P\i^y  proceed 

to  calculate  P22  and  p  i3=P3i  which  are  needed  as  inter¬ 
mediate  steps  in  a  perturbative  calculation  of  P33  [see,  for 
example,  chains  (6)  and  (7)]. 

To  second  order  in  the  fields,  it  follows  from  Eqs.  (21a) 
and  (25)  that  the  population  P22  can  be  written 

P22  =P22(  0 )  +P22(  S  )C  +  [P22(  S  )]•  C  . 


which,  when  substituted  into  Eq.  (21a),  may  be  combined 
with  Eqs.  (13d)  and  (22)-(26)  to  give 


P22(0)  =  - 


72  (r'i2)^+^^  (r'i2)^+^i 


lF(v)  , 


^  YT+W'lv)  2r,2 

P22(5)=— - ; - ^ ^  1  + - ^  .  (28b) 

(Yi2~''^4  )(7l2+'‘^l  72  +  'S 

Similarly,  to  second  order  in  the  fields,  it  follows  from 
Eqs.  (21e)  and  (25)  that  p^^  can  be  written  as 

p,3=p,3(A-KA')+p,3(A3-hA')e‘'''''^“^'’ 

-l-p,3(A-(-A'_)e^''''''^+*’'' 

+p,3(A  +-fA'_)e“'C*+'’’-'*“'*“^'’''  ,  (29) 

which,  combined  with  Eqs.  (21e)  and  (22)-(26),  yields 


P,3(A-I-A’)  = 


[r',3-/(A-l-A')](7'|2-'^> 


where 


(30a) 
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p,3(i  +  S'_)=- 


-xWW(y) 


[y',3-j(S++S')](y'2-/A+) 

-X*X-W{v) 

tr'u-'(S+A'-)](y',2-'A)  ’ 


(30b) 


Pl3' 


,(i++A'_)  = 


[y',3-/(A+  +  A'_)](y',2-'2^+) 


(30d) 


The  next  step  is  to  calculate  P23  =  P32-  To  third  order 
{30c)  in  the  fields,  it  follows  from  Eqs.  (2  Id),  (27),  and  (29)  that 
P23  can  be  written 


P23=P23(A')+p23(A'_)e-'<*'‘^  +  *'''  +  P23(A'-S)e-'''‘*-'‘''+p23(A'  +  g)e''*R-S'' 

+P23(  A  -B)e  +  +§)g.((«--r-)  R-(6  +  5'),l 

When  this  is  substituted  into  Eq.  (2 Id)  and  Eqs.  (27)  and  (29)  are  used,  one  may  obtain 
■<T'*P22(0)  +  JTPi3(^  +  ^  ')  +  'T+Pi3(A+  +A  ') 


P23(A')=- 

P23(A’_)  = 


P23(A'-S)  = 

P23(  A  '  +  S)  = 


723 -'A' 

-iY  *P22(0) + 'XP,3(  A  +  A  L ) + /AT +p  ,3(  A  + A '_ ) 
yi3-<AL 

-Y*p22(g)  +  «TPi3(A++A') 

r'2j~i(A'-S) 

-YV22(S)*+'T+P,3(A  +  A') 


P23(A'_-S)  = 
P23(AL+S)=- 


7'23  A  '  -hS) 

-Y-p22(^)  +  <TPl3(A.f+A'_) 

r23~<(A'_-S) 

-Y,V22(S)*  +  »7^P13(A  +  A'_) 
j'2}  —  i(A +  B) 


which  are  to  be  combined  with  Eqs.  (28)  and  (30). 

Finally,  we  calculate  P33.  To  fourth  order  in  the  fields, 
it  follows  from  Eqs.  (21b)  and  (31)  that  P33  can  be  written 
as 

P33  =  P33(0)+P33(S)e-'<*‘'-«'’ 

+P33(S')e''*'’^+«''’ 

+P33(B-S  ')e  -'((*+»''R-(»-S'id 

+P33(S  +  S')e"'''*'*'’‘^''®'"®'’'l  +  c.c.  (33) 

If  Eq.  (30)  is  substituted  into  Eq.  (21),  and  Eq.  (31)  is 
used,  one  can  obtain 

P33(0)=— -[/T'p23(A')  +  (;t-P23(A'_)]  ,  (34a) 

73 

p33(  5 )  = - )  /,v  P2 ,( A  '  -  8 )  -  ;,v'  *  [P2 ,( A  '  -I-  S )  ]  * 

73  +  '0 


P33(S')  = 


+  (T'  P23(A'_  -8) 
-,V[p23(A'_+8)]*)  , 

/y'-,P23(A')-iy'*[p23(A'-)]* 

73-1-/8' 


(34b) 


(34c) 


(31) 

(32a) 

(32b) 

(32c) 

(32d) 

(32e) 

(320 


(g  g,)-  YP23(A  --8)-Y,np23(^'  +  s')]* 
}'3+/(5— S') 


(34d) 


P33' 


(g  )  g.)-,  Y-P23(A'-8)-Y*[P23(A'-+S)]* 

r3+/(S+S') 


(34e) 


which  is  to  be  combined  with  Eqs.  (32),  (28),  and  (30). 

Equations  (25)-(34)  are  quite  general.  By  combining 
these  equations,  it  is  possible  to  show  that  P33(0)  has 
terms  which  vary  as  ItT'IM  V  f  IM7  +  7'  P.  'hat 

P33(8)  has  terms  which  vary  as  77*  It'I^  or  77+ 1 7'- I  ^ 
that  P33(8')  has  terms  which  vary  as  l7p7'-7'*  or 
I7 1 1 Y  7’*;  that  pyyib  —  b')  varies  as  77*  7'7'*,  and  that 
P33(8  +  5')  varies  as  77*7  *7’  •  In  this  paper  we  are 
concerned  only  with  the  contribution  to  p,,  that  varies  as 
the  product  of  the  four  distinct  field  amplitudes.  Experi¬ 
mentally,  this  contribution  could  be  isolated  by  field 
modulation  techniques.  The  contributions  to  p,,  which 
vary  as  the  product  of  the  four  amplitudes  are  represent¬ 
ed  by  the  terms  pyy{b~B')  and  p33(8  +  6')  in  Eq.  (34) 
(plus  their  conjugates)  [see  also  chains  (6)  and  (7)].  It  is 
possible  to  render  the  p3,(5  +  8')  contribution  negligible 
by  choosing  an  appropriate  phase-matching  condition.  If 
the  wave  vectors  k,k',k,,k'  ,  and  frequencies  8  and  5' 
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are  chosen  such  that 

ic=k+-k=-ic'=-(k'_-k')  , 

(35a) 

8=8’  , 

(35b) 

and  if,  in  addition. 

kLl  »1  or  |8tobsl  »1  , 

(36) 

where  L  is  the  cell  length  and  /gbs  **  some  observation 
time  in  the  experiment,  the  contribution  to  pjj  from  the 
_ I 


P33(5  +  5')  component,  which  varies  as  expj— /[(if 
— ir')-R  — (8  +  6')t])  [see  Eq.  (33)]  goes  to  zero  when 
averaged  over  the  cell  length  or  observation  time.  In 
what  follows,  we  assume  that  conditions  (35)  and  (36)  are 
valid,  and  we  calculate 

S(S)  =  p33(S-S)  +  [p33(S-S)]*  .  (37) 

From  Eqs.  (34d),  (32e),  (32d),  (28),  (30a),  (30d),  (35),  (23), 
and  (24),  one  finds 


5(5)  = 


XX%X'X-my) 

2r33+y3 


2r 


12 


1  +  - 


y2  +  /8 

+(r,3+r,3-r33) 


/[(yi3  +  /A'_)(yi3-/A')] 
1 


+  - 


(2r,3  +  y3) 

(y;3)^  +  (A  +  A')^ 


+  ■ 


1 


[y',3-/(A  +  A')](yi3-/A')  [y',3  +  /(A  +  A ')](y^3  +  /A '_ ) 


+  C.C. 


(38) 


r 


One  notes  immediately  that,  in  the  absence  of  collisions 
(r,y=0),  only  the  two-photon  resonance  centered  at 
A  +  A'=0  remains,  a  result  that  was  predicted  in  Sec.  II 
on  the  basis  of  conservation  of  energy  arguments.  When 
collisions  occur,  CARE-type  resonances  centered  at 
S.'=0,  A  1=0,  as  well  as  a  PIER4-type  resonance  cen¬ 
tered  at  8=0,  appear.  If  a  PIER4-type  resonance  exists, 
it  necessarily  follows  that  a  CARE-type  resonance  also 
exists.  This  interdependence  reflects  the  common  origin 
of  the  two  processes. 

To  examine  the  resonance  centered  at  8=0,  we  consid¬ 
er  the  limit  of  large  detunings  for  fields  E  and  E+,  such 
that 

|AUlA+UlA|^|A+!»/c«  .  (39) 

Moreover,  we  take  A'  =  0  in  order  to  maximize  the  signal, 
and  limit  8  to  values 

\8\'Sy2«ku,k'u  .  (40) 


KU  «y2  , 


(42) 


which  can  be  satisfied  for  nearly  parallel  beams  [the 
beams  must  be  slightly  misaligned  to  satisfy  inequality 
(36)].  When  inequality  (42)  holds,  one  finds 


<S(8)>  = 


XX*+X'X- 


X 


_ IPTv) _ 

[y23  +  /(kl-v  +  8)](y^3-/k'-v) 


X2r,2 


l-f 


y3+2r23 

y2  +  /8 


-l-c.c. 


(43) 


At  high  pressures  where  r23/y23 » 1,  for  |8|  <  y2  the  sig¬ 
nal  varies  as 


(5(8))a: 


xxXx'x'* 

A' 


ri2r23 
y2  +  /8 


+  C.C.  , 


(44) 


With  these  values  of  detunings,  the  second  and  third 
terms  of  Eq.  (38)  (and  their  conjugates)  can  be  neglected 
and  one  finds  a  velocity-averaged  signal  (denoted  by 
<■•■>), 


(5(8)>  = 


XX*+X'X'* 
y.,A'  ■' 


X 


_ IF(v) _ 

[y2,  +  /(kl  ■v-l-8)](y23  — /k’-v) 


that  is,  as  the  square  of  the  pressure.'^ 

An  important  feature  to  note  is  that  the  signal  varies  as 
TAffTyl*-  Thus,  for  a  nonvanishing  signal,  the  fields 
must  be  relatively  coherent.  This  result  is  in  contrast  to 
PIER4  where  the  signal  varies  as  the  absolute  square  of 
the  input  fields.*'  ^  Since  the  signal  calculated  above  de¬ 
pends  on  the  field  amplitudes,  it  is  possible  to  adjust  the 
relative  phases  of  the  fields  to  give  either  an  absorptive  or 
dispersive  shape  to  the  excitation  profile  (44). 


x2r  12 


1-1- 


y3+2r23 

y2  +  i(b-K-\) 


IV.  CALCULATION  OF  PIER  SIGNAL: 
DRESSED-ATOM  APPROACH 


-fc.c.  (41) 

To  avoid  any  significant  residual  Doppler  broadening 
of  the  8  =  0  resonance,  one  must  arrange  to  have 


We  now  outline  a  dressed-atom  approach  (DAP)  to  the 
problem  discussed  in  Sec.  III.  Each  radiation  field  is  now 
taken  to  be  quantized  rather  than  a  classical  field.  The 
Hamiltonian  for  the  system,  in  the  rotating-wave  or  reso- 
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nance  approximation  (i.e.,  |A|  «a);  |A'|  «(y')  is 

/f  =  Au<T 22  + )tr 33+  2  + 


7  =  1 

-1^2  .  (45) 

7  =  1 

where 

c7-V~‘^  +  ~|2)0|  ,  (46a) 

a^+=or^+  =  |3><2|  ,  (46b) 

aL=(<T\f  (46c) 


are  operators  that  act  on  atomic  states,  aj  and  Qj  are 
creation  and  annihilation  operators  for  field  mode  J,  and 
fj  is  a  coupling  constant  which  is  specified  below.  It  has 
been  assumed  that  each  field  consists  of  a  single  mode. 
The  fields  are  labeled  by  the  indices  1,2, 3, 4  which  corre¬ 
spond  to  the  fields  of  Fig.  4  as  follows: 


E,(k,,fti,ni )  — E(k,fl)  , 

(47a) 

E2(k2,ft2.”2)=E+(k+,fi  +  8)  , 

(47b) 

E3(k3,n3,n3)  =  E'(k',n') 

(47c) 

E4(k4,n4,«4)  =  E'_(k'_,n'-8')  , 

(47d) 

where  n,  is  the  average  number  of  photons  in  field  i. 

A  convenient  basis  set  consists  for  the  kets  is 
Iy;/j,,7j2;'i3<‘”4^  where  y  =  1,2,3  specifies  the  atomic 
state  and  n,  the  number  of  photons  in  field  i.  With  this 
set  of  basis  kets,  the  Hamiltonian  has  an  infinite  number 
of  nearly  degenerate  states  (e.g.,  all  states  with 
-fn2  “const,  are  nearly  degenerate).  However,  in  the 
limit  of  weak  fields,  the  number  of  nearly  degenerate 
states  which  are  of  importance  is  truncated.  If  our  goal  is 
to  calculate  P33  to  lowest  order,  the  relevant  basis  states 
and  associated  energies  are 

|l;n  )  =  |l;n,,n2'."3>”4)>  £i=0  (48a) 

|2;n  )  =  |2;n,  —  l,n2:”3.”4 ^2^^^  (48b) 

|3;n  )  =  |2;ni,n2  — I;n3,n4),  £'3=^(A  — 6)  (48c) 

I4;n)  =  |3;n,-l,n2;«3-l.«4>.  £4=«(A  +  A')  (48d) 

l5;n)  =  |3;n|  — 1,/J2;n3,n4  —  1),  £5=fi(A-l-A'-f5') 

(48e) 


|6;n)  =  |3;n|,n2~ll”3“l."4).  £6=^(A-I-A'  — 8) 

(481) 

l7;«>  =  |3;n,,M2-l:n3.n4-l>,  £7=^(A-l-A'  +  6'-8)  . 

(48g) 

An  interaction  representation  is  used  in  which  the  energy 
of  state  of  I  l;n  )  is  set  equal  to  zero  for  each  n.  The  label 
n  refers  to  a  set  of  specific  values  for  n ,  ,n2'”3>'’4-  Effects 
of  atomic  motion  are  neglected,  but  could  be  included 
easily  by  replacing  A,  A',  5,  and  8'  by  A,  A ',  S,  and  S'  of 
Sec.  II. 

Within  the  subspace  generated  by  the  basis  states 
I  1;m  ),...,  |7;rt  ),  the  Hamiltonian  takes  on  the  form 
shown  in  Table  I.  The  entries  are  shown  in  frequency 
rather  than  energy  units  and  the  coupling  constant 

(49) 

is  evaluated  at  the  average  number  of  photons  in  field  j, 
assumed  to  be  a  large  number.  The  coupling  constant  gj 
can  be  related  to  the  Rabi  frequencies  of  Sec.  Ill  as  fol¬ 
lows: 

g3  =  -/Y'e''‘ g*= -iX'-e‘'‘ '  ^  ■  (50) 

The  object  of  the  calculation  is  to  evaluate  the  upper 
state  population  defined  as 

P})=  2  (i\n,m;p,q\p\y,n,m;p,q) 

n,m,p,q 

=  2  i  {j-,n\p\j'\n’)  ,  (51) 

n.n'j.j'  =  * 

where  p=\tl;}{rl>\  is  the  density  operator  for  the  system, 
and  \j;n  )  is  given  in  Eqs.  (48d)-(48g).  The  second  sum 
in  Eq.  (51)  is  restricted  to  values  of  n  and  n'  such  that  the 
corresponding  occupation  numbers  of  states  \j',n)  and 
\j';n')  are  equal.  For  example,  if  j—j',  then  7i,=n'|, 
n2  =  n'i,  ny  =  n'y,  n^  =  n\-,  if  y  =4  and  y'  =  6,  then 
'’3~'*3>  't4='*4.  etc.  It  is  a 
straightforward  calculation  to  find  pyy  using  perturbation 
theory  for  the  Hamiltonian  of  Table  I.  That  calculation 
is,  in  effect,  identical  to  that  of  the  bare-atom  approach. 
In  this  section,  however,  we  wish  to  use  a  DAP  in  which 
a  part  of  the  Hamiltonian  is  diagonalized.  Using  weak- 
field  dressed  states  (i.e.,  dressed  states  correct  to  order  g. 


TABLE  I.  Subblock  of  matrix  H  (in  frequency  units)  in  the  |/ 

;ni.”2;n3 

,^4 )  basis. 

ID 

12) 

|3> 

|4) 

|5> 

|6> 

|7> 

<ll 

0 

ig* 

ig* 

0 

0 

0 

0 

<21 

-'gi 

0 

igt 

igt 

0 

0 

<31 

-igi 

0 

A-8 

0 

0 

ig* 

igt 

<4| 

0 

-ig} 

0 

A  +  A' 

0 

0 

0 

<51 

0 

-‘g4 

0 

0 

A  +  A'  +  6’ 

0 

0 

<61 

0 

0 

-ig} 

0 

0 

A-i-A'- 

6  0 

<7| 

0 

0 

-ig* 

0 

0 

0 

A-(-A’  +  8'-8 
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or^2)  defined  as'^ 

|/4;n>  =  |l;/i)+0Bl2;n)+0cl3;n>  ,  (52a) 

|B;n>  =  -0Sll;n)  +  |2;n)  ,  (52b) 

lC;n)=-6>J-|l,M)  +  |3,n)  ,  (52c) 

|4;n  >, l5;n  ), |6;n  ),  |7;n  )  as  in  Eq.  (48)  ,  (52d) 

where 


the  subscripts  n  and  n'  which  specify  the  quantum  num¬ 
bers  n^,n2,n■j,n^,  do  not  appear  in  matrix  elements  of 
It  is  assumed  implicitly  in  Eqs.  (55)  that  the  atom- 
field  interaction  is  turned  on  in  a  time  that  is  long  com¬ 
pared  with  I  A|  With  such  an  adiabatic  turn  on  of  the 
interaction,  terms  involving  time  derivatives  of  the  field 
amplitudes  are  neglected.  The  initial  values  for  density- 
matrix  elements  are  those  at  t  =  —  oo ,  before  any  atom- 
field  interaction.  Explicitly,  these  initial  values  are 


0B=ig,/A  , 


(53a) 


Pa  a  -p\ 


0c  =  '.?2/(^~S)  ,  (53b) 

we  replace  the  bare-state  Hamiltonian  of  Table  I  by  the 
corresponding  dressed-state  Hamiltonian  //^  shown  in 
Table  II.  The  Hamiltonian  of  Table  II  consists  of  two 
block-diagonal  sections  (separated  by  a  frequency  of  or¬ 
der  |A|)  which  are  coupled  by  the  atom-field  interaction. 
The  energy  levels  of  some  of  the  dressed  states  are  shown 
in  Fig.  5. 

In  the  absence  of  collisions,  the  time  evolution  of  the 
density  matrix  is  given  by 

p=— (54) 

where  the  tilde  indicates  the  interaction  representation. 
It  follows  from  Eq.  (54)  and  Table  II  that  density-matrix 
elements  evolve  as 


iH^ip,  -p  ,  Hfi)  , 

/=4 


p,v=  A,B,C  (55a) 


PLv=-i<^)vP!v-i  2  Hi^p  , 

nn 

H=  A.B.C 


v=A,B,C,  /=4,5,6,7,  (55b) 


LC  ' 


-KOirp, 


A,B,C 


where 


/,/'=4,5,6,7  ,  (55c) 


(56) 


It  is  to  be  noted  that,  neglecting  variations  of  gy  with  n. 


p,j^,^0,  i,j^A  .  (57) 

We  are  interested  in  a  perturbation  solution  to  Eqs. 
(55).  Since  is  of  order  g^,  and  Hjl  and  are  of  or¬ 
der  g ',  one  can  infer  from  Eqs.  (55)  that,  in  order  to  cal¬ 
culate  P33  to  fourth  order  in  g,  we  require 

p,v~0(g*)  , 

Pi^~0(g^)  , 

PiA-0{g^)  , 

(58) 

P^y-~0(g^)  , 

PAf,~Oig')  , 

PAA~0(g°)  , 

where,  unless  noted  otherwise,  we  adopt  the  convention 
that 


p,v  =  B,C  , 

(59) 

/,/'  =  4,5,6,7  . 

To  this  order  in  g,  Eqs.  (5)  become 


(60a) 

(60b) 

(60c) 

(60d) 

p,^,.--,co,rp,^,. 


TABLE  II.  Subblock  of  matrix 

(in  frequency  units)  in  the  dressed  basis. 

lA> 

\B) 

10 

14) 

|5> 

16) 

17) 

(-41 

0 

0 

0 

iOsgt 

iff*Bg: 

iff^g: 

<B\ 

0 

A 

0 

ig* 

ig: 

0 

0 

<cl 

0 

0 

A-6 

0 

0 

ig* 

ig: 

<4| 

-ig} 

0 

A  +  A' 

0 

0 

0 

(51 

~'g* 

0 

0 

A-)-A’-l-6’ 

0 

0 

<6| 

-iPcg} 

0 

-ig) 

0 

0 

A  +  A'-6 

0 

<71 

-iOcg* 

0 

-ig4 

0 

0 

0 

A -I- A'-)- 8' -6 
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with  initial  conditions 

PA„A^.=P°n‘  all  Other  py  =0  .  (60g) 

Relaxation  must  now  be  incorporated  into  Eqs.  (60). 
In  the  bare-atom  picture,  relaxation  is  accounted  for  sim¬ 
ply  by  writing 

( fj >rda*  =  - ( y )p fj  =  -y'jpfj  , 

r„=o,  r,=o.  (61) 

In  the  bare-atom  picture,  collisions  and  spontaneous 
emission  result  in  a  decay  of  the  various  density-matrix 
elements,  but  no  coupling  between  density-matrix  ele¬ 
ments  (collisions  are  not  energetic  enough  to  couple  the 
bare  states  and  the  repopulation  of  states  1  and  2  via 
spontaneous-emission  results  in  contributions  to  p^^  of  or¬ 
der  higher  than  g'*).  On  the  other  hand,  collisions  and 
spontaneous  emission  will  lead  to  coupling  of  dressed- 
state  density-matrix  elements.  To  arrive  at  the  relaxation 
equations  in  the  DAP,  (i)  Eqs.  (52)  are  used  to  write 
dressed  density-matrix  elements  p®  in  terms  of  the  bare 
ones  pfj-,  (ii)  Eq.  (61)  is  used  to  write  the  p*  in  terms  of 
the  pfj,  and  (iii)  the  inverse  of  Eqs.  (52)  is  used  to  reex¬ 
press  the  pfj  in  terms  of  the  p,y.  When  this  procedure  is 
carried  out,  the  resultant  expressions  are  very  complicat¬ 
ed.'^  However,  in  the  perturbation  theory  limit  of  Eq. 
(58),  they  take  the  relatively  simple  form 


r',}+iw,^ 


+  -I  X  — - ^ - 

.  _  •  _ " 

v=B,C  723  +  '";^ 


y23  +  ‘(% 


a—  A,B,C 


73  +  '"//' 


(63d) 

(63e) 

(63f) 


These  equations  are  exactly  equivalent  to  Eqs.  (34)  and 
lead  to  exactly  the  same  results  obtained  in  Eq.  (38). 

Although  Eqs.  (63)  are  “exact”  (in  the  perturbation 
theory  limit),  they  are  clearly  more  complicated  than  the 
corresponding  equations  in  the  BAP.  It  is  known  that 
DAP  expressions  simplify  significantly  in  the  so-called 
secular  approximation,'’"’  that  is,  when 


7^2«|A|  . 


(64) 


In  that  limit  all  the  terms  in  square  brackets  in  Eq.  (63) 
can  be  dropped  to  lowest  order  in  712/A,  and  one  obtains 
the  simplified  equations 


)rclax“ 

=0, 

(62a) 

0  , 

(65a) 

Op  a, 

—yi2PA„,„-o;y 

i2Pa„  > 

(62b) 

P..  »/  “ 

le.dirA 

(65b) 

v^./^^^relax 

2^^^jr'i2Pa,  A^~ 

+  (7'i2-72)(0^, 

,)  ,  (62c) 

iHFApln' 

7'i3  +  '"/a 

(65c) 

Op,^ 

Velax” 

■-r\2Pi„A„. . 

(62d) 

v=B,C  7  23  +  '<y/>i 

(65d) 

Op/„ 

^relax 

-7i3Pu,,  +  (7',3- 

-7  23)^MP/„a„. 

> 

HFp^  ,  -p,^^HFr 

(62e) 

_  n  H  "  " 

'a^A.BX  73  +  '<W/r 

(65e) 

(8p,  /a.' )„!«=- 73P, 


Equations  (60)  and  (62)  are  now  combined  and  the  fol¬ 
lowing  steady-state  solutions  are  obtained: 


P  A^A^.  P  nn 


;y0 

=  ^  (63b) 

7'i2  +  '"aM 

2e,e:r,jpl. 

72  +  '"mv 

+  [(7i2-72)(^IPm„ .  (63c) 


In  this  form,  the  role  of  collisions  is  readily  apparent.  In 
the  absence  of  collisions  (F,  =0),  Pbb~Pcc~Pbc 
~PcB~^  [Eq.  (65b)].  Conseque,.tly,  p,B=p,(-=0  from 
Eq.  (65d)  and  only  p/^ ,  corresponding  to  the  two-photon 
resonance,  contributes  lo  p,/  in  Eq.  (65e).  Thus,  in  the 
DAP,  the  absence  of  CARE-  and  PIER4-tyne  resonances 
can  be  attributed  to  the  fact  that  states  \B )  and  |C )  are 
not  populated  when  collisions  do  not  occur.  The  fact 
that  p^v~(^  (p,  v  =  fl,C)  in  the  absence  of  collisions  can  be 
connected  with  energy  conservation  —  there  is  no  mecha¬ 
nism  to  compensate  for  the  energy  defect  ^A  between 
states  I  A  )  and  |B )  (or  |C)).  Collisions  possess  frequen¬ 
cy  components  (of  order  of  the  inverse  collision  duration) 
to  drive  the  A~*B,C  transitions. 

As  a  specific  example,  we  e'aluate  Eqs.  (65)  for  the 
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same  set  of  parameters  used  in  arriving  at  Eq.  (38)  of  Sec. 
III.  We  set  8'=8  and  k=—k',  and  consider  only  the 
component  of  ,,  which  varies  as  XX\X'X'-  (or. 

equivalently,  as  gigJgjgJ  )  or  its  conjugate.  From  Eqs. 
(65)  and  Table  II,  it  follows  that  p^  ^  and  p^  4  have  this 

X  dependence.  In  the  interaction  representation  that  was 
chosen,  an  energy  subtracted  from  each 

_ I 


nearly  degenerate  subblock  of  H^.  This  implies  that  the 
total  density-matrix  elements  are  related  to  those  of  the 
interaction  representation  by 


Pi  I'M'^'^Pi  I'. 


-i  2 
;  =  i 


(66) 


Using  Eqs.  (66),  (48d),  (48g),  and  (47),  one  can  write 


P4,7,,(r)=<3;n,  -l,n2;M3-l,n4lp|3:n',,/i2~l;«3>"4-l ) 

Xexp(-/[(n,-w;)n-f(n2-n2)(n-f8)-l-(n3-n3)n'  +  (n4-n;)(n'-8')]fl  , 


(67) 


with  an  analogous  equation  for  P74.  In  calculating  the  contribution  to  P33  from  Eqs.  (67)  and  (51),  one  must  take  only 
those  terms  which  are  diagonal  in  the  photon  numbers,  that  is,  terms  for  which 


«'i=n|  —  1,  n2=n2  +  (>  n\=n^  +  l 


(68) 


For  these  values  of  n  and  n',  the  phase  factor  in  Eq.  (67)  vanishes  identically  when  8  =  8'.  Consequently,  the  contribu¬ 
tion  to  P33  from  P47  and  p-,^,  denoted  by  5(8)  as  in  Sec.  Ill,  is  given  by 


S(8)  = 


(p4  78.  .8  .  ,8  .  8  ,  )-t-c.c. 


(69) 


"r"2’"3'"4 

Using  Eqs.  (65)  and  (50)  and  Table  II,  one  can  obtain 


5(8)  = 


xxlx'x- 

y^A{A~8) 


2r 


12 


1  + 


y3+2r23 

yj  +  iS 


-f  ■ 


r3+2r,3 


[r^3  +  j(A'-8)](ri3-/A')  (y',3)^-)-(A-l-A')^ 


(p°>+c.c.  , 


(70) 


where 

(p°)—  X  (fl)  lp”|n|  -  1  )<n2lp2i"2  + ( ^ 

X  <«3ip3|n3  —  1  )<n4lp4|n4-)- 1 )  (71) 

and  p”  is  the  free-field  density  matrix  for  field  1.  ” 

There  are  two  points  we  wish  to  note  concerning  Eq. 

(70) .  First,  we  see  that  it  reproduces  all  but  the  last  term 
of  Eq.  (38).  The  last  term  of  Eq.  (38)  corresponds  to  a 
nonsecular  term  in  the  DAP  which  is  not  included  in  the 
secular  approximation  to  the  DAP,  Eqs.  (65)  [this  term  is 
included  the  full  DAP  equations  (63)].  The  nonsecular 
term  is  negligible  if  one  has  (A-fA')~0  or  A'~0,  but 
would  be  comparable  to  the  other  terms  if  both 
lA-f  A'l  »ku  am-  :A'|  »k'u.  Second,  it  is  interesting  to 
note  the  dependence  on  photon  statistics  given  in  Eq. 

(71) .  For  coherent  states  .f  each  fi'*’-!,  the  average  densi¬ 

ty  matrix  <p°)  is  equa'  ...  nity.  However,  if  even  one  of 
the  fields  has  a  total!.-  r  j  .n  phase  (i.e.,  one  field  is  in  a 
pure  number  state;,  <  .  anishes.  As  mentioned  in  Sec. 

Ill,  the  fields  must  b.  relati-  .,  coherent  for  a  nonvanish¬ 
ing  5(8). 

Thus,  from  eitht  a  bare-  or  dressed-atom  approach, 
we  arrive  at  the  same  result — only  the  two-photon  reso¬ 
nance  is  present  in  the  absence  of  collisions.  The  DAP 
offers  a  unified  picture  of  the  pressure-induced 
resonances — CARE-type  resonances  result  from 
collision-induced  population  of  dressed  states  \B  )  and 
|C),  while  PIER4-type  resonances,  centered  at  8  =  0,  re¬ 
sult  from  collision-induced  creation  of  the  pgc  coherence. 


I 

There  remains  one  fly  in  the  ointment,  however.  Popu¬ 
lations  Pbb  and  pec-  well  as  coherence  pBc<  are  not  ex¬ 
actly  zero  in  the  absence  of  collisions  owing  to  contribu¬ 
tions  from  nonsecular  terms  arising  from  spontaneous 
emission.  These  nonsecular  terms  would  lead  to  CARE- 
or  PIER4-type  resonances,  even  in  the  absence  of  col¬ 
lisions,  terms  which  are  known  to  vanish  identically  from 
Sec.  III.  Of  course,  if  one  solved  the  full  set  of  DAP 
equations  (63),  he  would  find  that  these  terms  vanish  as  a 
result  of  interference  of  various  pathways  in  the  perturba¬ 
tion  calculation.  However,  the  simple  interpretation  in 
terms  of  the  DAP  would  be  tarnished  somewhat  by  these 
nonsecular  terms.  The  situation  can  be  remedied  if  com¬ 
plex  dressed  states  are  introduced  (see  Appendix).  In  that 
case,  one  can  show  that  pgg,  p^c,  and  p^B  are  identically 
zero  in  the  absence  of  collisions — the  simple  interpreta¬ 
tion  of  the  DAP  is  retained. 

V.  SUMMARY 

By  studying  a  three-level  atom  interacting  with  four  ra¬ 
diation  fields,  we  have  shown  how  both  collisionally  aided 
radiative  excitation  (CARE)  and  pressure-induced  extra 
resonances  in  four-wave  mixing  (P1ER4)  can  be  explained 
in  a  unified  manner.  Using  both  a  bare-atom  picture 
(BAP)  and  dressed-atom  picture  (DAP)  of  the  reaction, 
we  were  able  to  argue  that  both  types  of  pressure-induced 
extra  resonances  (PIER)  could  be  explained  on  the  basis 
of  conservation  of  energy.  In  particular,  we  have  found 
the  DAP  to  be  a  convenient  vehicle  for  viewing  these 
pressure-induced  effects. 
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In  the  near  future,  we  plan  to  extend  this  work  in  two 
directions.  First,  using  both  the  BAP  and  a  new,  semi- 
classical  dressed-atom  approach,  we  will  calculate  PIER 
signals  arising  in  (i)  fluorescence  beats,  (ii)  a  “four-level” 
atom  interacting  with  four  radiation  fields,  and  (iii^  two- 
photon  ionization.  Second,  we  will  use  both  the  quan¬ 
tized  and  semiclassical  DAP  to  reanalyze  pressure- 
induced  effects  produced  in  four-wave  mixing  on  a  “two- 
level”  atom. 


</<1/4  0  =  (B"|E")  =  1  , 

(A^\B^)  =  (B\A^)=0.  (A5) 

To  lowest  order  in  the  fields,  the  eigenkets  and  eigenfre- 
quencies  of  H  and  are 

M0  =  |1>+(Pb!2>,  |B0  =  |2>-^5I1)  (Abat 
|7‘'>  =  |T>+^fli2>,  |50  =  |2)-(pS|T>  (A6b) 
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APPENDIX:  COMPLEX  DRESSED  STATES 


—0,  Wg  =a)g*  =( A  — j/yj)  (A7) 

where 


<Pb  ~ 


‘8\ 

A-i/ra 


>  fa  ~~ 


‘g\ 

A +  1/72 


(A8) 


and  |l  ),|2)  and  |l  ),|2)  are  eigenkets  of  H  and  re¬ 
spectively,  in  the  absence  of  the  applied  field.  Note  that 
<Pa~^B  [Eq-  (53a)]  when  y2-~0- 

An  arbitrary  ket  It/i)  is  expanded  (in  the  interaction 
representation)  as 


As  discussed  in  Sec.  IV,  it  would  be  more  satisfactory  if 
dressed-state  populations  pgg  and  pcc  vanished  identical¬ 
ly  in  the  absence  of  collisions,  not  only  in  the  secular  ap¬ 
proximation.  It  is  possible  to  arrange  for  this  to  occur  if 
“complex”  dressed  states  are  used.  To  illustrate  the  con¬ 
cepts,  we  consider  a  two-level  atom  interacting  with  a 
single  quantized  field  (e.g.,  field  1  of  Sec.  IV).  In  the  ab¬ 
sence  of  collisions,  the  state  amplitudes  and  a  2  for 
states  1  l;n  )  and  |2;n  —  1 ),  respectively,  evolve  as 

(Ala) 


<d2  =  ~'(y2/2)a2  +  Aa2“'giai  .  (Alb) 


where  the  energy  of  state  |l;n)  has  been  set  equal  to 
zero.  These  equations  imply  that  a  complex  Hamiltonian 
(in  frequency  units) 


0  /gf 

-/g,  A-i/y2 


(A2) 


can  be  introduced  for  which 


2 


idj=  y  =  1,2  . 

1  =  1 


(A3) 


We  now  define  complex  dressed  states  (CDS’s)  as  eigen¬ 
states  of  H.  In  the  spirit  of  the  perturbative  calculation 
of  this  paper,  the  CDS’s  are  r  btained  to  first  order  in  g , , 
but,  in  principle,  the  CDS’s  can  be  defined  for  arbitrary 
g,.  Since  H  is  non-Hermitian,  some  care  must  be  exer¬ 
cised.'* 

Let  1^4“^)  and  be  eigenkets  of  H.  Since^/f  isjiot 
Hermitian,  {  A’^\B'^)=^0,  in  general.  Define  \A 
as  eigenkets  of  H  ,  having  eigenvalues  which  are  the 
complex  conjugates  of  those  of  states  \A^)  and  In 

other  words. 


HW)  =o)a\a.‘^)  , 

//^|aO=<|aO,  a=.4,B.  (A4) 

The  eigenkets  of  H  and  form  an  orthonormal  set'* 


\i^)=a^^\A^)+a%\B‘^}  ,  (A9) 

where 

=  a=  A,B  .  (AlO) 

The  density  matrix  (in  the  interaction  representation)  is 
defined  as 

p  =  |^><^|  ,  (All) 

with  density  matrix  elements  in  the  complex  basis 
defined  by 

.  (A12) 

The  time  evolution  of  density  matrix  elements  is  obtained 

from  p  =  (/B)“'(f/p— pf/^).  Taking  matrix  elements  as 
in  (A  12),  one  finds 


^aP=-'-2:[2/5,plp-p^x(2/g,)*]  ,  (A13) 


where  the  completeness  relation,  =  1,  has 

been  used.  Defining  a  complex  dressed  Hamiltonian 
as  one  having  matrix  elements  between  states  (a ‘‘I  and 
la*^),  one  can  use  Eqs.  (A2)  and  (A6)  to  write 


H. 


D  — 


0 

0 


(A14) 


From  Eqs.  (A  13)  and  (A  14),  it  immediately  follows  that 


p'^AA^^  ,  (A15a) 

P’^AB~P  BA~  ^  ^^P'^AB  >  (A15b) 

P%B~  ~yiP'^BB  ■'  (A  15c) 

subject  to  the  initial  conditions 

P  AA~^’  P  AB~P  AB~P  BB~^  ■  (A16) 


In  the  absence  of  collisions,  but  including  spontaneous 
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emission,  Pbb~P^ab~Pba~^  this  complex  dressed- 
state  basis.  It  remains  to  include  collisional  relaxation  in 
Eqs.  (A  15).  The  procedure  for  carrying  this  out  is  out¬ 
lined  in  Sec.  IV.  In  terms  of  the  complex  dressed  states, 
one  finds,  to  the  order  indicated  in  Eq.  (58), 


P'^AA  • 

(A  17a) 

P‘AB~~(\y'l~^P'>P'^AB~^\l‘P*BP’^AA  > 

(A  17b) 

PbB~  ~y2PBB'^^n^<PBP‘^AB'^‘PBpBA  ^ 

(A  17c) 

These  equations  can  now  be  generalized  to  the  three- 
level  problem  of  the  main  text.  The  Hamiltonian 
of  Table  II  can  be  replaced  by  in  which  the  6^  of 

Eq.  (53)  are  replaced  by 

<Pa=  \  't-  .  P  =  B,C  (A18a) 


P'AA=P°n 


(A  19a) 


'  W  V  . 
“rt  n 


r\2 


Afi 

xO 


(A19b) 


■  +  ^l2^<P*P'f^„A+(P^^A„ 


- '  2  - 

yi3'^‘^IA  ii  =  B,C  y  n'^ ‘^lA 


— ;  V' 

f^FAP\,„. 

^  2, 

v=B,C 

y23  +  '"/;4 

-i  S 

^F'p" - 

— 

a^A.B.C 

(A  19c) 
(A19d) 


(A19e) 


(A19f) 


and  the  6*  of  Eq.  (53)  are  replaced  by 

,  P^B,C  ■  (A18b) 

Using  Eq.  (A  13),  one  can  show  that  the  steady-state  solu¬ 
tions  analogous  to  Eqs.  (63),  correct  to  the  order  of  Eqs. 
(56),  are 


where  is  still  defined  by  Eq.  (56).  In  the  absence  of 
collisions  (Ty  and  p^Jp,v  =  B,C)  are  identical¬ 

ly  zero.  Consequently,  p]b~P‘^ic~^^  the  sole  contri¬ 
bution  to  P33  arises  from  the  two-quantum  resonances. 
With  collisions  present,  both  the  stepwise  and  two- 
quantum  resonances  contribute.  In  the  absence  of  col¬ 
lisions,  the  vanishing  of  p^,.  (p,v  =  BC)  is  no  longer  re¬ 
stricted  to  the  secular  approximation.’’ 
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to  first  order  in  g,  and  g2.  A  term  —  6*|l;n|  —  l,n2 
-l-l;n,,n4>  should  be  added  to  Eq.  (52b)  and  a  term 
—  0JI  l;n,  -(-  l,n2  ~  I:u3,n4 )  to  Eq.  (52c).  It  then  follows  that 
|B;ni,n2;n,,n4)  =  |C;ni -I,n2  +  l:”3.'i4)  and  that  all  the 
dressed  states  can  be  labeled  simply  by  1 /I ,n2;n  1,04 ), 
Ifl;n',,n2;”3.n4 ).  for  arbitrary  which  satisfy 

ni+n2=n  \  -f  n  j.  The  equations  of  the  text  are  unaltered  if 
these  "correct"  dressed  states  are  used.  The  physical  argu- 
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ments  given  in  this  paper  are  based  implicitly  on  a  transition 
from  an  initial  density-matrix  element 

to  a  final-state  density-matrix  element 

<3;n,  —  l,n2;n3~l.n4lpl3:nl,n2“Un3■n4“l)• 

Such  a  final-state  density-matrix  element  does  not  contribute 
to  the  final  state  population  p^.  To  be  rigorously  correct,  the 
physical  arguments  should  be  based  on  a  transition  from  an 
initial-state  coherence 

0;ni-l-l,/i2;”3  +  l.”4lpll;'>i.'i2  +  l;n3.n4+l) 
to  a  final-state  population 
<3;n,,/i2:«3.«4lp|3;ni,n2:«3.«4)  • 


Mathematically,  the  two  approaches  can  be  mapped  into  one 
another.  We  have  chosen  the  first  method  to  explain  the  re¬ 
sults  since  the  physical  arguments  can  be  presented  in  a  some¬ 
what  simpler  fashion  using  that  method. 

'■^P.  R.  Berman  and  R.  Salomaa,  Phys.  Rev.  A  25,  2667  (1982). 

•’The  Jif’s  appearing  in  Eq.  (70)  have  no  n  dependence  since  the 
g's  have  been  defined  at  n  =  n  [Eq.  (49)]  on  the  assumption 
that  h  is  large  for  each  mode.  The  n  dependence  of  the  ^"’5 
could  be  included  by  using  Eqs.  (49)  and  (50)  with  h  set  back 
equal  to  n  in  Eq.  (49). 

'*P.  M.  Morse  and  H.  Feshbach,  Methods  of  Theoretical  Physics 
(McGraw-Hill,  New  York,  1953),  pp.  884-886. 

•^Complex  dressed  states  have  also  been  used  to  calculate  the 
collision-induced  Rayleigh  gain  for  an  ensemble  of  two-level 
atoms;  H.  Friedman  and  A.  D.  Wilson-Gordon,  in  Advances 
in  Laser  Science  III,  edited  by  A.  C.  Tam,  J.  L.  Gole,  and  W. 
C.  Stwalley  (AlP,  New  York,  1988),  pp.  121-129. 


